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Processes such as double Drell-Yan and same-sign WW production have contributions from double 
parton scattering, which are not well-defined because of a 5^ (z_i_ = 0) singularity that is generated 
by QCD evolution. We study the single and double parton contributions to these processes, and 
show how to handle the singularity using factorization and operator renormalization. 



I. INTRODUCTION. 

Single parton scattering (SPS) processes such as Drell- 
Yan, pi +P2 — > l + 1~ , shown in Fig.JT] involve one parton 
in each hadron colliding via a hard interaction. Factor- 
ization allows one to write the cross section as 



da DY 



CO 



AX\AX2 Xi X2 



[f q {xi)h{x 2 )+f^ Xl )f q {x 2 )], (1) 



at leading order. Here, <To = <iira 2 Q q /(3Q 2 ) is the short- 
distance qq —> partonic cross-section, Q is the 
lepton-pair invariant mass which sets the hard scale, and 
Q q is the quark charge. The f(xi) are the (single) par- 
ton distribution functions (PDFs), and the momentum 
fractions fixed by the invariant mass and rapidity 
of the lepton pair. The transverse momenta of the lep- 
tons is integrated over. If the transverse momenta of the 
leptons is measured, the expression in Eq. ([I]) is modi- 
fied and involves transverse-momentum-dependent PDFs 
instead. 

















? \ 




V 



FIG. 1. Leading order diagram for single Drell-Yan. The 
cross-section is shown as the cut of a forward-scattering am- 
plitude. The incoming hadrons are double lines, the gauge 
bosons are the wiggly lines, and the final state leptons are 
the cut loop. The hard interaction is given by shrinking the 
gauge boson lines to a point. 
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FIG. 2. Double parton scattering contribution to double 
Drell-Yan. The two hard interactions, given by shrinking the 
gauge boson lines to a point, are vertically separated in the 
figure by z±, which is not constrained by measurements. 



One can also have processes where two partons in one 
hadron collide with two partons in the other hadron, 
which is known as double parton scattering (DPS). The 
leading-order DPS contribution to the double Drell-Yan 
cross section is shown in Fig. [2] and is (schematically) 
given by [1] 



..DPS 



dxida^da^da^ 



<7q J d 2 zj_F(xi,x 2 ,zj_)F(x 3 ,x 4 ,z J _) 



(2) 



The measured lepton-pair invariant masses and rapidities 
fix the momentum fractions Xi, and the transverse mo- 
menta are again integrated over. The hard double-parton 
scattering cross section er 2 is of order (7ra 2 /Q 2 ) 2 , i.e. the 
same order as the square of (To in Eq. (HI) . Precise numer- 
ical factors can be found in Refs. [21 [3]. F(xi,Xj, z±) are 
the double PDFs (dPDFs), which depend on the mo- 
mentum fractions Xi,Xj and the transverse separation 
zj_ ~ 1/Aqcd of the two hard collisions. The transverse 
separation Zj_, or its Fourier-space analog k^, is not de- 
termined by the measurement and must be integrated 
over m Eq. ^, even if one measures the transverse mo- 
menta of all the leptons. The dPDFs F(xi, Xj,z±) are of 
order Aq CD , and the z j_ integral of two dPDFs is also of 
order Aq CD , so the DPS cross section is of order Aq Cd ctq 
and Aq CD /Q 2 suppressed relative to SPS. 

There are additional terms that contribute to Eq. ^ 
with spin and color correlations [5J [5] and interference ef- 
fects [3 E] , which involve soft functions. Explicit expres- 
sions for their contribution can be found in e.g. Ref. [3], 
and are briefly discussed in the appendix. The color cor- 
relation and interference contributions are Sudakov sup- 
pressed, and thus negligible at high energies [7] . 



II. MIXING OF SINGLE AND DOUBLE PDFS. 

The intuitive description of DPS is somewhat mis- 
leading because there is also a contribution from SPS 
to double Drell-Yan, which mixes with DPS under the 
renormalization group evolution. For example, the graph 
shown in Fig. [3] leads to a mixing of the gluon PDF f g 
with the qq dPDF F q q. This contributes to the renormal- 
ization group evolution of the color-summed dPDF F^g 



2 




FIG. 3. Mixing of the gluon single PDF f g with the qq double 
PDF F a a. 



as [HIEJEE] 

dFUx 1 ,x 2 ,z_ L ) 



d// 



7T 



P qg (x) = T F [x 2 + {l-xf] 



fg{x 1 +X 2 ) 
Xl+X 2 



(3) 



in terms of the usual g — > q splitting function P qg (x). 
The ". . ." denote non-mixing contributions given in 
Refs. [3, 5, 6, 8, 9], and are not important for the present 
discussion. The S^ 2 \zj_) form for the mixing arises be- 
cause the gluon splits into a qq pair at the same point. 

We now discuss the modifications to Eq. ([3| for the 
spin correlated, color correlated and interference dou- 
ble PDFs, as defined in Ref. [3]. The mixing anoma- 
lous dimension for F AqA - is the same as for F qq , and for 
FsqSq ^ is given by replacing x 2 + (1 — x) 2 in Eq. ([3| 
by — 4x(l — x). For the interference double PDFs I qq 
and I AqAq the x 2 + (1 — x) 2 is replaced by — 2x(l — x), 
and for I} q5q it is replaced by 2[x 2 + (1 — x) 2 ]. For the 
color-correlated dPDFs F T and I T the color factor T F is 
replaced by {C F ~ C A /2)T F . 

As was pointed out in Refs. [5J [5J, the <5^ 2 )(zj_) form 
of the mixing contribution in Eq. ([3| is problematic. 
Under RG evolution F^(xi, x 2 , Zj_) develops a <5^ 2 ^(zj_) 
contribution, or equivalently, its _L-Fourier transform 
F^(xi,X2,kj_) develops a kj_-independent contribution. 
The integral of this term with the other dPDF in Eq. ([2| 
gives 

y"d 2 z_L (5 (2) (z_L)/ 3 (a;i + :r 2 ) -^(2:3, x 4 , z_i_) 



= —fg( x i + x 2 )F q 1 Jx 3 ,X4,z 1 _ = 0) 



(4) 



However, there is also the convolution of the mixing con- 
tribution to both dPDFs with each other 

dV f g (xi + x 2 )5 ( - 2 \z x )f g (x 3 + x 4 )6^(z ± ) 

= f g (x 1 +x 2 )f g {x 3 +x i )5( 2 \z ± =0)=7, (5) 

which is singular. In momentum space, this singularity 
arises from the d 2 kj_ integral of a constant. In this letter 
we show how this problem is solved by a careful QCD 
factorization analysis of a physical cross-section, such as 
double Drell-Yan. One hint is provided by the observa- 
tion that the <5*- 2 )(zj_) contribution is when the qq pair 
are at the same z± , which overlaps with the single parton 
scattering region. 
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FIG. 4. Gluon contribution to deep inelastic scattering. The 
dashed area is the hard interaction, and is shrunk to a point. 
The wiggly lines are photons, and the springy lines are gluons. 




FIG. 5. Single parton scattering contribution to double Drell- 
Yan. 



III. DIAGRAMMATIC ANALYSIS. 

The Drell-Yan cross-section in Eq. (JlJ follows from fac- 
torizing the QCD graph for the cross-section, shown in 
Fig. [T] as the imaginary part of the forward scattering 
amplitude. To study (and define) DPS we look at a 
physical process, such as double Drell-Yan. This receives 
both SPS and DPS contributions. Due to the mixing 
in Eq. ([3]), the separation of the cross section into SPS 
and DPS contributions depends on the renormalization 
scheme and renormalization scale /1. This is analogous to 
deep-inelastic scattering, where the cross-section can be 
written in terms of quark and gluon distributions. The 
total is well-defined, but the split depends on renormal- 
ization (or factorization) scheme and scale \i. The same 
graph in Fig. [3] gives the mixing of the quark and gluon 
PDFs (left figure) or the one-loop jg hard-scattering 
cross section (right figure). These individual pieces are 
scheme dependent. 

We now discuss the various contributions to the dou- 
ble Drell-Yan cross-section, starting with the single par- 
ton scattering contribution in Fig. [5] The cross-section 
is [a/(47r)] 2 suppressed relative to Eq. 0, but is still 
leading twist. The leading double parton contribution is 
shown in Fig. [2] and leads to Eq. The mixing contri- 
bution is shown in Fig.[6j This graph has a short-distance 
contribution to double Drell-Yan of the form 



a ~ 7^0 fg( x l + x 2 )F^ Jx 3 ,x 4 ,z_ L =0). 



47T 



(6) 



where the entire diagram is shrunk to a point. This same 
diagram also has a contribution that corresponds to the 
mixing graph in Fig. [3] [in analogy with Fig. [2j. From 
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FIG. 6. Mixing contribution to double Drell-Yan. 





FIG. 9. The [/ d 2 z ± F(z x )F(z ± )] dPDF operator. 
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FIG. 10. The [f g F± q (z± = 0)] operator. 



FIG. 7. Chirality suppressed mixing contribution to double 
Drell-Yan. 



Eq. Q we see that it is of the same size as Eq. Q. Fi- 
nally, there is also the mixing graph shown in Fig. [7J 
which leads to two twist-three PDFs, which are qqg ma- 
trix elements. However, these are chirality suppressed by 
light quark masses, and can therefore be neglected. 

The problematic term in Eq. Q arises from the dou- 
ble mixing graph in Fig. [8j This graph is down by 
[a s /(47r)] 2 [a/(47r)] 2 compared to single Drell-Yan, and 
is leading order in the twist expansion. It is of the same 

order in coupling constants as the [S^ (z_i_)l term that 
arises from mixing. The part of the diagram that would 
produce Eq. |5]) is a quadratically divergent kj^ integral. 
The two powers of k_L convert the leading- twist f g f g 
gluon PDFs into the Aq CD /Q 2 -suppressed dPDF contri- 
bution, since |kj_| ~ Aqcd- However, this bound on kj^ 
does not appear in perturbative calculations, and there 
is no physical scale of order Aqcd that enters the graph. 
The quadratically divergent scaleless integral vanishes in 
dimensional regularization, so the contribution in Eq. ^ 
is absent and there is no problem! 
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FIG. 8. Double mixing contribution to double Drell-Yan. 



IV. OPERATOR RENORMALIZATION. 

The above discussion can be rephrased in terms of the 
renormalization of composite operators. The dPDF is 



[F(z ± )] , 



(7) 



keeping only the _L argument explicit. The square brack- 
ets emphasize that we are using the renormalizcd four- 
quark operator (defined explicitly in Refs. OEJE]) to get 
the finite dPDF. Fig. [5] gives a contribution 



Jd 2 z±F(* ± )F(* ± ) 
where the key point is that 



(8) 



^ / dV [F(z ± )} [F(z ± )} . (9) 



d 2 z ± F(z ± )F(z ± ) 



The integral over zj_ includes the point zj_ =0, and one 
needs to perform an additional renormalization for prod- 
ucts of operators at the same space-time point. This 
is analogous to the well-known result that [</> 2 (a;)] ^ 
[4>{x)\ [4>{x)\ in an interacting theory. 

We represent the dPDF operator in Eq. ([£]) by Fig. [9] 
The mixing graph in Fig. ([6| gives 



d/i 



J d 2 z ± F* ? {x-i , x 2 , z j_ ) F± q (x 3 , x A , z ± ) 

IT "\Xi+X 2 J 



(10) 



fg(x 1 +X 2 ) X , ' 

— ; Fg (x 3 , x A , z ± = 

X\ + x 2 



p„ 



x 3 



19 \ 

7T VX3+X4 



FgJx 1 ,x 2 ,z ± = 0) 



fg{x 3 +X i ) 



X 3 + X4 

where the right-hand side is the matrix element of the 



operator shown in Fig. 10 Both sides of the equation 
are operator matrix elements of order Aq CD . 
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FIG. 11. Diagram for mixing the [/ s i^ 5 (zi = 0)] operator 
with gluon PDFs f g f g . 




The mixing of \f g F^ a (z± = 0)] into two single gluon 
PDFs shown in Fig. [IT] vanishes, 



^[/,J&(»l = 0)]=0. 



(11) 



The left-hand side is order Aq CD , whereas the single 
gluon PDFs are order Aq CD . There is no order Aqcd 
_L observable to compensate for the dimensions, so the 
two objects cannot mix. There is no longer any prob- 
lem with a S^(zj_ = 0) as in Eq. Explicitly, the 
computation gives the integral 



d 2 " 2e k ± 



(12) 



This is the same integral that arises in Eq. QSJ) . The dif- 
ference is that Eq. ( 12 ) is evaluated in fractional dimen- 



sion using dimensional regularization, where it vanishes, 
whereas Eq. ^ arises after renormalization, and is eval- 
uated in integer dimension, where it is singular. The non- 
mixing contribution for the color-correlated dPDF F T is 
also modified at z± — 0, as discussed in the appendix. 

The above discussion shows that the singular quan- 
tity in Eq. ([5| does not enter the renormalization group 
evolution of the dPDF. A recent paper by Gaunt and 
Stirling |10j discusses the double parton scattering singu- 
larity in gg — ¥ ZZ, which is an equivalent way to study 
Fig. [8] They suggest that one should drop the f g f g mix- 
ing term in Eq. |5]), in agreement with our result. A 
similar conclusion was also reached in Ref. [llj . 



V. ESTIMATING THE SIZE OF VARIOUS 
CONTRIBUTIONS. 

We thus have a well-defined factorization formula for a 
physical process such as double Drell-Yan. The factorized 
cross-section has the following (schematic) form 



Clf q fg + ^f g fg + £ 3 lfgF(z± = 0)1+^3 [F(Z A 

f d 2 z ± F(z x )F(z x ) 



+C4 



0)/ 9 ] 
(13) 



where Cj are the partonic cross sections for the hard scat- 
tering, and contains both single and double parton scat- 
tering contributions. There are many possible terms in 
Eq. (131 for all the different possible diparton combina- 
tions, such as uu, ud, ug, etc. in the dPDFs. For dou- 
ble Drell-Yan, the ci_4 terms contribute [a/(47r)] 2 ao, 



FIG. 12. Leading single PDF contribution to W W pro- 
duction. 



[a/(4ir)) 2 [a s /(4ir)] 2 cr , [a s /(4 7 r)]A 2 JCD cr 2 and A 



respectively to the cross section, where the coupling con- 
stants are evaluated at the hard scale Q. 

Same-sign WW production is a process in which the 
single PDF contribution is suppressed. The dPDF con- 
tribution is the same size as for double Drell-Yan, from 
Fig. [2] However, the single PDF contribution requires 
two additional partons (jets) in the final state, as shown 
in Fig. [12] The c\ contribution to the cross section is 
~ [a/(47rsin 2 W r )] 2 [ Q; s/(47r)] 2 cro, where <j 1S now the 
single-W partonic cross section, leading to the naive es- 
timate ~ 10~ 9 ao- H only a factor of ir is included with 
each a, one finds instead ~ [a/(nsm 2 9w)] 2 [a s /^] 2 cro ~ 
1Q- J a . The dPDF contribution F^F^ to W+W+ pro- 
duction is of order Aq CD (t 2 . This process has been stud- 
ied in more detail in Refs. [T2HT5] . 

The evolution in Eq. ( 10 1 vanishes in this case be- 



cause ud cannot mix with a gluon. The cross section also 
gets contributions from other flavor combinations. For 
example, the gluon single PDF f g f g term contributes 
at order [a/(i% sin 2 #w)] 2 [a s /(47r)] 4 <7o and the F U uF^ 
dPDF term contributes at order [a s /(47r)] 4 AQ CD (7 2 . The 
<5( 2 )(OjJ problem is four higher orders in a s compared 
with double Drell-Yan. This mixing contribution is a 2 
suppressed relative to the f g f g contribution. 

In conclusion, we have studied the mixing between 
single and double PDFs, and shown how to solve the 
S i - 2 \0±) problem. The factorization theorem for a phys- 
ical process such as double Drell-Yan has both single 
and double parton contributions. These mix, and so are 
only separately defined after a choice of renormalization 
scheme. DPS is not a physical process by itself, but one 
contribution to a physical process. 
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Appendix A 

For the color-correlated (and interference) contribu- 
tions, there is also a subtlety at zj^ = for the non- 
mixing contributions. As discussed in Refs. [3J OH], the 
effect of soft radiation is nontrivial in these cases, and is 
described by a soft function depending on zj_ ~ 1 / Aqcd • 
As zj_ — > 0, the soft radiation becomes energetic and is 
part of the short-distance physics. The soft function dis- 
appears, since S(z± = 0) = 1. This is reflected in a 
change in the anomalous dimensions at zj_ =0, which 
we illustrate below for the color-correlated dPDF F^-. 

At zj_ = 0, the Wilson-line diagrams I A' and IB' in 
Ref. [3j are changed to I A' — —IB' = I A, yielding 



/./,- 



dF£ j (x 1 ,x 2 ,0±_) _ a, 

d[l 7T 



dz 

z 



P w{— )F qq (z,x 2 ,0 x ) 



Pj«{x) = Cf(^^) + - \c a (\ - x) . (Al) 

The first term in Pj q is the usual q —> q splitting function, 
the second term causes the dilution of color correlations 



d*j£(*) = -% 



(A2) 



The evolution in Eq. (JAlJ) is quite different from the ex- 
pressions for zj_ =/= in Ref. [3,. In particular, for z± =/= 
the dPDFs and soft functions have rapidity divergences, 
which leads to an additional evolution involving a new 
renormalization scale v [TBI [T7] . 
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